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Abstract — Multiphoton  nonlinear  effects  and  Raman  processes  of  Rydberg  systems  (Rydberg  atom  gas  or 
Rydberg  matter)  by  intense  laser  fields  are  considered.  Some  experiments  and  basic  elements  of  quantum  theory 
of  Rydberg  matter  are  briefly  reviewed.  A correlation  between  experimental  data  and  theoretical  description  is 
made. 


1.  INTRODUCTION 

Systems  consisting  of  Rydberg  atoms  (RAs)  attract 
much  attention  during  the  last  decade.  On  the  one  hand, 
both  a single  isolated  Rydberg  atom  and  RA  clusters 
are  of  considerable  importance  for  laser  physics  and 
technologies.  On  the  other  hand,  modem  laser  physics 
and  technologies  provide  new  opportunities  for  the 
experimental  investigation  of  many  problems  that  are 
of  fundamental  significance  for  quantum  physics  and 
quantum  optics.  Some  of  such  studies  have  been 
already  performed.  Among  those,  we  should  mention 
investigations  around  the  borderline  between  classical 
and  quantum  physics  [1],  analysis  of  information  con- 
tent and  transfer  of  quantum  states  [2,  3],  and  experi- 
ments with  single  atoms  and  ions  (see  [4]  and  refer- 
ences therein). 

Creation  of  a coherent  superposition  of  Rydberg 
states  proved  to  be  very  useful  for  the  investigation  of 
the  classical  limit  of  quantum  mechanics.  The  spectrum 
of  Rydberg  levels  is  nearly  equidistant  for  large  values 
of  the  principal  quantum  number.  In  this  sense,  an  RA 
forms  the  same  coherent  states  as  a quantum  oscillator. 
An  RA  in  a strong  laser  wave  or  a strong  magnetic  field 
is  an  attractive  object  for  studying  quantum  systems 
whose  classical  analog  exhibits  a chaotic  behavior 
[5,  6].  Nevertheless,  at  the  moment,  there  is  no  satisfac- 
tory description  of  the  correspondence  between  the 
quasi-classical  and  purely  classical  behavior  of  an  RA. 
However,  an  attempt  was  made  recently  to  describe  the 
motion  of  wave  packets  consisting  of  Rydberg  states 
with  very  large  orbital  momenta  (circular  Rydberg 
wave  packets)  within  the  framework  of  the  so-called  de 
Broglie-Bohm  pilot  wave  theory  [7]. 

Modem  experimental  methods  make  it  possible  to 
produce  RA  systems  with  sufficiently  high  concentra- 
tions, when  the  treatment  of  such  a system  as  a gas 
becomes  incorrect.  Rydberg  matter  (RM)  is  usually 
understood  as  a system  of  Rydberg  atoms  where  the 


mean  distance  between  the  atoms  is  comparable  with 
RA  sizes.  Due  to  the  collectivization  of  motion  of  RA 
electrons  [8],  RM  is  reminiscent,  to  some  extent,  of  the 
metal  phase  of  matter  with  an  ultralow  gas  density  [9]. 
Rydberg  matter  is  a metastable  condensed  phase  of  the 
nonequilibrium  state  of  a medium  resulting  from  the 
decay  of  a nonideal  low-temperature  plasma  or  arising 
under  the  laser  excitation  of  atomic,  molecular,  cluster, 
or  aerosol  gas  species  to  Rydberg  states.  The  review 
[10]  provides  the  description  of  physical  and  chemical 
properties,  as  well  as  the  main  relaxation  processes  of 
RM  decay.  Currently,  there  are  several  well-established 
methods  of  RM  creation.  Experimental  studies  of  RM 
started  ten  years  ago  [11]. 

The  most  obvious  and  straightforward  technique 
implies  a two-photon  excitation  of  alkaline  and  alka- 
line-earth atoms  to  Rydberg  states.  Vitrant  et  al.  [12] 
made  the  first  attempt  to  ensure  conditions  meeting  the 
requirements  of  Mott  phase  transition  in  an  RA  gas. 
Rydberg  atoms  with  principal  quantum  numbers  n = 
30-45  were  produced  in  a laser-pumped  cesium  vapor. 
In  accordance  with  the  Mott  criterion,  the  phase  transi- 
tion to  the  metallized  state  was  expected  for  RA  con- 
centrations corresponding  to  the  close  packing  of  RAs, 
i.e.,  for  RA  concentrations  meeting  the  inequality 
NCs(n2aBr )3  > 1,  where  NCs  is  the  concentration  of  Cs 
atoms  excited  to  the  nth  Rydberg  state.  Hence,  we  see 
that  the  critical  concentration  NCs  scales  as  rr6,  while 
more  accurate  estimates  predicted  that  the  transition  to 
the  metal  phase  should  be  observed  even  with  the  ratio 
of  the  mean  interatomic  distance  to  the  atom  size  equal 
to  three.  The  experimentally  measured  dependence 
scaled  as  nr4.  Subsequent  studies  have  shown  that,  for 
mean  distances  between  cesium  RAs  no  less  than 
15  times  the  RA  size,  an  RA  gas  rapidly  decays  in  an 
avalanche-like  fashion  through  a long  series  of  sequen- 
tial chemical  reactions  involving  the  formation  of  vari- 
ous neutral  and  charged  molecular  complexes.  There- 
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fore,  the  required  concentrations  were  not  achieved 
with  the  use  of  this  method,  and  no  RM  was  produced. 
The  authors  of  [1 2]  argued  that  RM  could  be  created  by 
considerably  reducing  the  influence  of  relative  RA 
motion,  which  could  be  done,  for  example,  by  substan- 
tially lowering  the  temperature  (below  7=1  K)  or  by 
drastically  decreasing  time  intervals  of  RA  generation 
and  detection  (down  to  the  picosecond  scale). 

It  was  found  out,  however,  that  mixtures  of  different 
sorts  of  RAs  (the  so-called  RM  alloys)  can  be  produced 
in  a rather  simple  way  by  using  thermal-emission 
diodes.  These  methods  were  thoroughly  investigated  in 
a series  of  studies  by  Holmlid  and  his  colleagues  (see 
[13]  and  references  therein).  This  group  of  researchers 
developed  a method  allowing  the  generation  of  high 
concentrations  of  RAs  using  special  graphite  mem- 
branes as  cathodes  in  thermal-emission  diodes.  Cesium 
gas  passing  through  the  pores  of  a membrane  heated 
with  IR  laser  radiation  up  to  temperatures  T = 1450°C 
was  transformed  into  a mixture  of  a low-temperature 
nonideal  plasma  and  a gas  consisting  of  different  sorts 
of  RAs.  This  approach  allowed  the  necessary  RA  con- 
centrations to  be  achieved,  permitting  the  conditions 
required  for  the  condensation  of  RAs  into  RM  to  be 
provided. 

The  general  theory  of  condensed  excited  states 
(CES),  which  was  developed  earlier  in  [9],  was  applied 
to  a specific  gas  of  cesium  RAs  [14],  and  the  theory 
describing  the  decay  of  RM  through  different  relax- 
ation channels  (involving  spontaneous  emission  and 
Auger  processes)  was  developed  [15, 16]. 

An  “optical  lattice”  consisting  of  cold  RAs  trapped 
in  a periodic  potential  [17]  is  another  interesting  RA 
system,  which  is  intensely  studied  at  the  moment. 
Atoms  in  an  optical  lattice  possess  an  evident  flexibility 
and  can  be  employed  for  the  generation  of  quantum 
states  and  subsequent  control  of  quantum  states  [18],  as 
well  as  for  quantum  computations  [19].  Such  RA  sys- 
tems substantially  increase  their  efficiency  due  to  the 
large  number  of  Rydberg  states  and,  which  is  even 
more  important,  due  to  large  decoherence  times  attain- 
able in  RAs. 

We  should  note  also  that  a considerable  progress  has 
been  achieved  recently  in  the  generation  of  strongly 
cooled  neutral  plasmas,  including  ultralow-temperature 
plasmas  [20-23]  and  plasmas  allowing  RA  formation 
[24-26], 

We  should  specially  mention  unique  experiments  on 
the  formation  of  an  ultracold  plasma  [23,  26].  The 
decay  of  such  a plasma  is  one  of  the  possible  ways  to 
produce  RM.  An  ultracold  neutral  plasma  with  the 
mean  electron  density  on  the  order  of  ne  = 1 09  cm'3  was 
created  through  the  photoionization  of  Xe  atoms  with 
the  density  na  = 1010  cm-3  cooled  with  laser  radiation 
down  to  ultralow  temperatures  (Te  = 100  mK,  Ta  = 
10  pK).  Three-body  recombination  is  suppressed  with 


such  plasma  parameters,  and  the  kinetic  (thermal) 
energy  of  plasma  species  is  lower  than  the  mean  Cou- 
lomb energy,  which  makes  the  plasma  a strongly  Cou- 
lomb-coupled system.  Xenon  atoms  were  excited  to  a 
metastable  state  (with  a lifetime  of  43  s)  for  laser  cool- 
ing through  the  transition  with  X,  = 882  nm  with  subse- 
quent trapping  in  a magneto-optical  trap  with  the  size 
d ~ 200  pm.  About  106  atoms  were  trapped  in  this  trap. 
Two-photon  ionization  at  ^ = 514  nm  (the  green  beam) 
was  employed  for  this  purpose.  The  difference  A 
between  the  final  level  of  photoelectron  excitation  and 
the  ionization  potential  of  Xe  was  varied  within  the 
range  of  A/k  = 0.1-1000  K by  changing  the  frequency 
of  the  green  beam.  The  lower  limit  was  imposed  by  the 
bandwidth  of  this  beam  (0.007  cm-1).  Thus,  the  initial 
velocity  of  a photoelectron  was  v = [2A Jm)m.  The 
number  of  photoions  was  controlled  by  changing  the 
energy  of  laser  pulses.  The  maximum  value  of  this 
energy  was  1 mJ  per  10-ns  pulse,  which  allowed  up  to 
2 x 105  ions  to  be  produced.  The  Debye  screening 
length,  / = {kT/4ne2ne}U2,  could  reach  the  minimum 
value  l ~ 500  nm,  which  implies  that  the  condition  of 
plasma  formation  is  satisfied,  i.e.,  I < d.  The  opposite 
inequality,  l > d,  is  characteristic  of  the  ionized-gas 
phase,  rather  than  the  plasma  phase.  The  equality  / = d 
determines  the  critical  concentration  of  the  phase  tran- 
sition. Killian  et  al.  [23]  argued  that  the  energy  of  a 
trapped  electron  in  the  gas  phase  is  completely  trans- 
formed into  the  kinetic  energy  in  the  plasma,  with  Te  <, 
A/k,  when  the  system  passes  through  the  critical  point, 
and  the  above-specified  values  of  Te  and  Ta  are  achieved 
as  a result  of  the  thermal ization  process.  It  can  be  easily 
seen  that  the  conditions  of  RM  formation  are  realized 
when  the  inequality  A < 0 is  met.  In  this  case,  a dense 
gas  of  highly  excited  supercooled  RAs  characterized  by 
a considerable  overlapping  of  electron  shells  may  be 
produced.  Consequently,  a Mott  transition  to  a metal- 
lized or  another  plasma-like  state  may  occur  under 
much  more  favorable  conditions  than  the  conditions 
attained  in  the  above-described  situation.  Many-body 
quantum  interactions  become  important  in  this  regime, 
prevailing  over  elastic  and  inelastic  collisions  in  gas 
phases.  Killian  et  al.  [26]  studied  RA  formation  in  an 
expanding  plasma  with  a temperature  varying  within 
the  range  of  1-1000  K and  densities  of  1 05— 10*°  cm-3. 

It  was  found  that  up  to  20%  of  the  initial  charges 
recombine  within  approximately  100  ps.  This  rather 
long  lifetime  is  inconsistent  with  the  mechanism  of 
three-body  recombination.  The  energy  of  RA  coupling 
is  approximately  equal  to  the  kinetic  energy  of  free 
electrons  surviving  through  the  recombination.  Thus, 
the  conditions  implemented  in  these  experiments  not 
only  allow  the  formation  of  a large  number  of  RAs  in 
an  ultracold  neutral  plasma,  but  also  permit  RM  forma- 
tion. 

Summarizing  the  aforesaid,  we  should  note  that,  on 
the  one  hand,  RM  exhibits  many  properties  characteris- 
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tic  of  metals.  On  the  other  hand,  several  unusual  prop- 
erties are  observed  in  this  case.  For  example,  the  den- 
sity of  RM  corresponds  to  rarefied  gases,  RM  is  trans- 
parent in  the  optical  range,  possesses  a high  plasticity 
and  a sufficiently  high  coefficient  of  surface  tension, 
which  allows  the  formation  of  a sharp  boundary 
between  an  RM  and  gas  media. 

In  this  paper,  we  provide  theoretical  estimates  for 
the  nonlinear-optical  properties  of  RAs  and  RA  sys- 
tems, including  RM,  responsible  for  multiphoton  pro- 
cesses. We  will  demonstrate  that  the  linear  and  nonlin- 
ear susceptibilities  for  RAs  and  RM  are  rather  high  as 
compared  with  typical  linear  and  nonlinear  susceptibil- 
ities of  ground-state  atoms.  This  conclusion  is  con- 
firmed by  the  results  of  experiments  devoted  to  the 
observation  of  Raman  scattering,  four-wave  mixing, 
and  multiphoton  photoionization.  We  will  then  provide 
a brief  review  of  some  of  the  RM  properties.  Theoreti- 
cal results  will  be  compared  with  the  results  of  recent 
experiments. 

2.  NONLINEAR  SUSCEPTIBILITIES 
OF  A RAREFIED  GAS  OF  RYDBERG  ATOMS 

We  will  illustrate  our  method  of  estimating  nonlin- 
ear susceptibilities  and  specify  our  main  approximation 
by  estimating  linear  susceptibilities. 

The  well-known  formula  for  the  linear  susceptibil- 
ity Xnl  (to)  of  an  arbitrary  quantum-mechanical  system 
residing  in  some  stationary  state  |/>  can  be  represented 
as 


,(a)  ,(P) 

®im  & mi 

hca  + hm„ 


.<P)  ,(«) 

aim  ami 

h(t>-hc o„ 


(1) 


where  d\ and  d^fj  are  the  matrix  elements  of  the 
operator  of  the  dipole  moment  vector  of  the  system 
under  study;  a and  P are  the  indices  of  Cartesian  coor- 
dinates x,  y,  and  z;  ha>  is  the  energy  of  the  laser  wave 
photon  with  a frequency  to;  ha>mi  = em  - e(  ; and  em  is  the 
energy  spectrum  of  the  system.  Summation  in  Eq.  (1)  is 
performed  over  the  intermediate  stationary  states  | m)  of 
the  system  with  m i.  The  factor  N for  an  arbitrary  sys- 
tem should  be  omitted.  In  the  case  of  a rarefied  gas,  the 
system  is  reduced  to  an  atom,  and  the  factor  N,  appear- 
ing in  front  of  the  sum,  is  understood  as  the  concentra- 
tion of  atoms  in  the  gas.  Expression  (1)  holds  true  also 
for  a crystal  medium  if  only  bound  electrons  are  taken 
into  consideration.  In  the  latter  case,  the  factor  N stands 
for  the  number  of  unit  cells  per  unit  volume. 

To  simplify  the  final  formulas,  we  assume  that  the 
electric  fields  of  laser  waves  are  linearly  polarized 
along  the  direction  of  the  x-axis.  Then,  it  is  sufficient  to 
restrict  our  consideration  to  a single  tensor  element 

Xxx  (co).  We  will  denote  this  quantity  as  x**  (to)  = %fl), 
omitting  the  coordinate  indices. 


In  what  follows,  Rydberg  states  will  be  denoted  as 
|n).  We  set  |i)  = |n)  and  represent  %(1)  as  a sum  of  two 

terms:  x(1)  = XrV  + xlll  • The  term  Xr,  is  given  by 
Eq.  (1)  with  |i)  = |n),  | m)  = |nj),  n*nx.  Thus,  we  need  to 
perform  summation  in  Eq.  (1)  only  over  the  Rydberg 
states  different  from  fixed  |«>.  In  contrast  to  the  first 

term,  the  second  term  Xrest  involves  summation  over 
m * nx.  In  other  words,  summation  has  to  be  made  in 
this  case  over  all  the  states  except  for  the  Rydberg  ones. 
We  will  consider  nonresonant  cases,  when  (D  (0ml  for 
all  the  discrete  states  m and  l. 

To  estimate  xi's! , it  is  sufficient  to  set  da  = ea  ( a is 
the  atomic  scale  on  the  order  of  the  Bohr  radius).  This 
approximation  holds  true  for,  at  least,  several  atomic 
transitions  with  maximum  values  of  dml.  Under  nonres- 
onant conditions,  all  the  denominators  of  the  form 
he o ± h(Om,  can  be  replaced  by  hma,  where  ma  is  some 
typical  atomic  transition  frequency  (h(Oa  = e2la  = 2Ry, 
where  Ry  is  the  Rydberg  constant).  Using  these 

approximations,  we  reduce  Eq.  (1)  to  Xrest  = 
-N{da)2{h(n^rl  = -Na3.  The  same  result  can  be 
obtained  if  we  estimate  the  susceptibility  of  a gas  with 

atoms  residing  in  the  ground  state,  x«°  = -Na3,  since 
|j)  in  Eq.  (1)  in  this  case  is  related  to  the  ground  state, 
and  summation  over  m = nx  gives  a negligible  contribu- 
tion. Thus,  we  finally  have  Xrest  = Xo°  • f°r  example, 
a ~ 10~8  cm  and  N ~ 1019  cm-3,  then  x^  ~ 10“5. 

Now,  we  use  a similar  approach  to  estimate  Xr/y  • 
The  electron  density  distribution  in  an  RA  reaches  its 
maximum  at  the  distance  on  the  order  of  rRy  ~ n2a. 
Therefore,  we  can  employ  the  following  order-of-mag- 
nitude  estimate:  dnn]  ~ n2ea  = dRy  (e.g.,  dnl_X  ni  = 

ea^n  Jti 2 - 12  ~ n2ea  for  Kn  [27]).  Then,  for  Rydberg 

1 e2  Ry  „ ...  . ^ 

levels,  e„  = r — = — ^ , we  find  that  to,,  „ _ x = — < 

2 n a n n 

to.  Substituting  co„,„_  i and  dn,u  ~ dRy  for  Xr^  and  dnn[ 
in  Eq.  (1),  we  arrive  at  the  following  estimate  for  co,,,,, 
linearly  depending  on  the  principal  quantum  number  n: 

.2 

XRy  = o„,„_,  = -2n{Na)  = 2 . (2) 

nco 

Hence,  we  find  that,  with  n > 1 , the  absolute  value  of 
X^y  is  much  greater  than  the  linear  susceptibility  XiP 

of  a conventional  gas  and  xist  = ■ Therefore,  we 

can  assume  with  a satisfactory  accuracy  that  the  linear 
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susceptibility  of  a Rydberg  gas  is  x(1)  = x^y  + Xrest  ~ 

XRy  • In  other  words,  we  can  completely  neglect  the 
contribution  of  all  the  intermediate  states  except  for  the 
Rydberg  ones. 

Expanding  Eq.  (1)  as  a series  in  a small  parameter 
to„iM  /go,  we  derive  a more  accurate  expression  for  xn): 


indices  and  arguments  will  be  omitted  below  for  the 
sake  of  convenience).  The  well-known  expression  for 
X(3)  is  written  as  [28,  29] 


,(3)  _ N 
1 h\_ 


p -0)dlm.d, 


imf*m , i 


xff  - x0>  = 


C0„ 


nn  ,1 


(3) 


2 N 
h($A 


iwv 


X { [(3(0  - O>(I1],)(2C0  - (0„,2,)((0  - to(„3,.)]' 

- K®  + tom1/)(2o)  - co„,2,)(co  - to,n3,.)]-'  (6) 

+ [(®  + co„,it)(2co  + to,„2,)(<»-  coWj,)]_1 


We  can  now  use  the  sum  rule,  which  is  well  known  in 

2 

atomic  physics  [27];  • Invoking 

this  rule  and  restricting  our  analysis  to  the  first  term  in 
Eq.  (3),  we  arrive  at  the  following  result: 


(i)  (i)  e2N 

ARy  “A.  — 

fflO) 


(4) 


This  expression  also  gives  the  linear  susceptibility  of  a 
rarefied  gas  of  free  electrons  in  the  classical  model. 
This  is  not  surprising,  since  the  coupling  of  an  electron 
with  an  ion  in  an  RA  is  weak  as  compared  with  the  pho- 
ton energy  in  the  optical  range.  For  the  convenience  of 
our  further  analysis,  we  will  represent  Eq.  (4)  in 
another  form.  We  introduce  the  electromagnetic  radius 
2 

£ 

of  an  electron  re  = — -2  and  the  wavelength  of  optical 
me 

field  A = X/2n  = c/co  to  arrive  at 


XRy  = X(l)  = -(A^A2).  (5) 

The  identity  of  Eqs.  (2)  and  (5)  implies  that  n = 
0.5reA2a~?l.  Substituting  the  values  re  = 10~13,  a = 0.5  x 
10~8  cm,  and  A = 10~5  cm  into  this  expression,  we  find 
that  n = 40.  This  result  clearly  shows  that  our  approxi- 
mations and  Eq.  (2)  for  Xr/  are  quite  adequate  for  reli- 
able estimates  on  n within  the  interval  of  n ~ 10-100. 

To  demonstrate  our  method  of  estimation  of  nonlin- 
ear susceptibilities,  we  will  consider  susceptibilities 

(3) 

XapyS  (®i>  ®2>  <%)  responsible  for  nonlinear-optical 
effects  of  the  third  order,  such  as  third-harmonic  gener- 
ation (THG),  four-wave  mixing,  self-focusing,  self- 
phase modulation,  two-photon  absorption,  stimulated 
Raman  scattering,  etc.  To  make  these  estimates,  we  will 
employ  the  approximations  substantiated  above  for  the 

estimation  of  xip  (<*>)•  In  particular,  to  describe  THG  in 
linear  polarized  laser  waves  with  the  same  polariza- 
tions,  we  should  estimate  %xxxx  (go,  go,  go)  = x(3)  (tensor 


- [(®  + tom,/)(2(o  + wm2,)(3co  + ay,,)]"1 }. 

This  susceptibility  can  be  also  represented  as  a sum  of 

two  terms:  x(3)  = XRy  + Xrest  • Each  of  these  terms  is 
given  by  Eq.  (6)  with  different  types  of  summation:  for 

XRy  * \mj)  - l«/)>  where  j = 1 , 2,  3 and  tij  is  the  principal 

quantum  number  of  Rydberg  states,  and  for  xrest . \mj)  * 
Using  the  same  assumptions  as  those  employed 

above  to  estimate  xl'i  . we  find  that  the  expression  for 

Xrest  coincides  with  the  relevant  nonlinear  susceptibil- 

ity  xfl  for  a rarefied  gas  of  ground-state  atoms: 


Xrest  = X?  = -~L  = -{Nay)%  = -Na3±  (7) 

fi V e Ea 


In  the  final  formula  (7),  Ea  = e/a 2 is  a typical  value  of 
the  intraatomic  electric  field.  Setting  Ea=  107  esu,  N = 

1019  cm-3,  and  a = 10  8 cm,  we  find  that  Xrest  = xT  = 

X(a\EaTx  = 10",2esu. 

Finally,  we  can  estimate  Xr^  using  the  same 

approach  as  that  employed  above  to  find  x^y  • Perform- 
ing simple  algebraic  calculations  involving  series 
expansions  of  expressions  of  the  form  of  Eq.  (6)  in  a 
small  parameter,  we  derive 


(3)  _ 1 NdR yC0„  „_i 

XRy  - ( 


V -3 


CO 

, 3,1  8 dg  1 

y (n C0j  n 


(0 


(8) 


, »t  3.  n 
” ~(Na 


where  we  used  the  notations  co  = co0,  dRy  = n2da,  and 
co„  „ _ , = co fl/n3  for  the  final  formulas. 
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(3) 

As  can  be  seen  from  the  final  expression  (8),  %Ry  = 

5 

%(3)  = X «3>^  > i-e->  Xity  *s  a power  function  (of  the  fifth 
power).  With  n = 10,  we  have  %R3)  ~ 10 4%^3)  = 10'8  esu. 

Thus,  xr3)  is  much  greater  than  %(3) . This  implies 
that,  for  example,  the  efficiency  of  third-harmonic  gen- 
eration in  a gas  of  RAs  with  a concentration  N = 
1014  cm-3  is  nearly  equivalent  to  the  THG  efficiency  in 
a gas  of  ground-state  atoms  with  the  concentration  Na  = 
1019  cm-3. 

Thus,  RAs  in  the  gas  phase  possess  very  high  non- 
linear-optical  susceptibilities,  which  is  confirmed  by 
the  available  experimental  data.  However,  we  should 
expect  that  RAs  may  exist  not  only  in  the  form  of  a gas, 
but,  although  it  may  seem  paradoxical,  in  a liquid  or 
even  a solid-state  phase  of  matter,  where  nonlinear  sus- 
ceptibilities are  also  very  high. 

3.  CONDENSED  PHASE  OF  RYDBERG  ATOMS 
(RYDBERG  MATTER) 

At  first  glance,  formation  of  Rydberg  matter  is 
impossible,  since  the  lifetime  of  an  excited  atom  is  very 
short — about  10~8  s.  Then,  such  an  atom  is  transferred 
to  the  ground  state  emitting  a photon.  An  excited  atom 
may  relax  even  faster  through  an  Auger  process  when  it 
collides  other  excited  atoms.  In  this  case,  an  electron  of 
the  first  atom  transfers  its  energy  to  an  electron  of 
another  atom,  which  becomes  ionized  as  a result  of  this 
collision,  while  the  first  atom  is  transferred  to  the 
ground  state. 

Nevertheless,  condensed  excited  states  may  arise 
under  certain  conditions,  and  the  existence  of  such 
states  is  currently  a well-documented  fact. 

The  idea  of  producing  a metallized  phase  in  a gas 
was  put  forward  back  in  1980  [30].  The  theory  of  CES, 
based  on  the  method  of  density  functional  [9]  (see  also 
references  in  [13]  and  review  papers  [10,  31]),  was 
developed  several  years  later. 

Formation  of  CES  and  RM  starts  with  the  condensa- 
tion process  caused  by  the  polarization  (van  der  Waals) 
mutual  attraction  of  Rydberg  atoms  toward  each  other. 
For  Rydberg  atoms,  this  attraction  is  n7  times  stronger 
than  for  nonexcited  atoms  or  molecules.  When  the  elec- 
tron shells  of  RAs  overlap,  i.e.,  for  interatomic  dis- 
tances comparable  with  RA  sizes,  exchange  interac- 
tions come  into  play.  These  interactions  determine  elas- 
tic repulsive  forces  and  the  mean  stable  distance 
between  RAs  in  metastable  RM. 

When  the  electron  shells  of  RAs  overlap,  an  elec- 
tron may  move  with  equal  probabilities  around  its  par- 
ent ion  and  around  the  neighboring  ions  before  it 
relaxes  to  a lower  lying  state.  Such  valence  electrons  in 
a cluster  are  usually  called  delocalized,  or  collectivized 


electrons.  The  area  of  motion  of  electrons  substantially 
expands  under  these  conditions.  The  sizes  of  this  area 
increase  as  quantities  multiple  of  the  number  of  atoms 
in  a cluster.  This  leads  to  an  even  sharper  decrease  in 
the  kinetic  energy  of  electrons  in  clusters  consisting  of 
RAs  relative  to  an  H2  hydrogen  molecule.  As  a conse- 
quence, the  coupling  energy  becomes  rather  large:  an 
RA  cluster  becomes  stable  with  respect  to  evaporation 
of  Rydberg  atoms.  These  effects  are  similar  to  the  for- 
mation of  metals  of  alkaline  and  alkaline-earth  atoms. 
In  this  sense,  a cluster  can  be  considered  as  a structure 
consisting  of  a metallized  condensed  phase,  i.e.,  the 
phase  of  RM. 

General  description  of  the  physical  properties  of 
RM  is  usually  provided  in  terms  of  a simple  model  of 
an  RA  cluster.  A cluster  in  this  case  is  understood  as  an 
assembly  of  RAs  closely  packed  into  a face-centered 
cubic  or  closely  packed  hexagonal  structure,  character- 
istic of  solid-state  crystals  or  the  near  order  in  liquids. 
An  elementary  model  consists  of  identical  RAs  resid- 
ing in  spherically  symmetric  S states.  With  large  quan- 
tum numbers  n,  spherically  symmetric  S states  still 
exist.  An  electron  cloud  with  n > 1 in  the  S state  radi- 
cally changes  in  this  case  from  the  electron  cloud  with 
n = 1.  In  the  considered  case,  the  electron  cloud  is  con- 
centrated within  a narrow  layer  with  a thickness  naRr 
near  the  surface  of  a sphere  with  the  radius  n2aRr.  An 
RA  in  the  S state  can  be  considered  as  a hollow  atom, 
similar  to  a bubble,  since  the  density  of  the  electron 
cloud  on  the  nucleus  and  in  most  of  the  inner  part  of  the 
atomic  volume  is  negligibly  small.  For  substantially 
asymmetric  states,  toruslike  RAs  are  usually  consid- 
ered. In  both  cases,  a minimum-size  cluster  consists  of 
seven  RAs  and  has  a shape  of  a planar  hexagon.  A min- 
imum three-dimensional  cluster  of  spherical  RAs  con- 
sists of  13  closely  packed  RAs.  All  the  valence  elec- 
trons in  RM  are  collectivized  and  are  concentrated  in 
thin  layers  along  the  boundaries  of  centered  unit  cells 
(Wigner-Seitz  cells).  Clusters  of  larger  sizes  are  con- 
structed by  adjoining  additional  atoms  to  minimum- 
size  clusters.  Spatially  extended  crystal  samples  with 
macroscopic  sizes  may  arise  in  this  manner. 

An  Ashcroft  pseudopotential  can  be  employed 
instead  of  the  Coulomb  potential  to  describe  the  inter- 
action between  an  electron  and  a nucleus  (or  an  elec- 
tron and  an  ion  core)  within  the  framework  of  the  con- 
sidered model.  The  Ashcroft  pseudopotential  elimi- 
nates the  singularity  characteristic  of  the  Coulomb 
potential,  replacing  this  singularity  by  a step  in  the  cen- 
tral area  of  an  RA,  where  the  nucleus  with  electrons  of 
inner  atomic  shells  is  located,  i.e.,  in  the  region  of  the 
ion  core.  In  the  first-order  approximation,  this  replace- 
ment does  not  have  a significant  influence  on  the  accu- 
racy of  calculations.  The  electron  cloud  of  a valence 
electron  has  a very  low  density  around  the  ion  core,  i.e., 
around  the  step.  Parameters  of  this  step  are  chosen  in 
such  a way  as  to  give  correct  values  of  energy  and  the 
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electron  density  distribution  in  an  isolated  RA.  The 
introduction  of  a pseudopotential  and  the  use  of  the 
method  of  density  functional  allows  a variation  proce- 
dure to  be  employed  for  the  determination  of  the  mini- 
mum energy  density  in  RM. 

This  theory  has  much  in  common  with  the  general 
method  of  calculation  of  physical  parameters  for  sim- 
ple metals.  The  difference  is  that,  in  the  case  of  simple 
metals,  e.g.,  alkaline  metals,  the  lower  energy  states  of 
the  ion  core  are  filled  with  electrons,  while  the  lower 
states  in  RM  are  vacant.  However,  this  difference  has 
no  influence  on  the  main  physical  and  chemical  proper- 
ties of  metals  and  RM,  since  these  properties  are  deter- 
mined by  valence  electrons.  The  presence  or  absence  of 
electrons  in  lower  levels  (inside  ion  cores)  influences 
only  the  lifetime  of  metastable  RM. 

The  CES  theory  allows  many  parameters  and  prop- 
erties of  RM  to  be  calculated  in  a rather  simple  way  or 
reliably  estimated  for  different  initial  principal  quan- 
tum numbers  n.  In  particular,  this  approach  can  be  used 
to  assess  interatomic  distances,  the  equilibrium  density, 
the  coupling  energy,  elasticity  modulus,  the  coefficient 
of  surface  tension,  the  speed  of  sound,  the  melting 
point,  the  boundary  of  the  transparency  range,  the 
Fermi  level,  electric  conductivity,  and  many  other 
parameters. 

The  most  important  achievement  of  the  above- 
developed  approach  relative  to  the  method  employed 
for  the  calculation  of  metal  parameters  is  associated 
with  the  fact  that  the  distribution  of  the  electron  density 
(e-density)  in  RM  is  quite  nonuniform.  In  simple  met- 
als, the  e-density  remains  virtually  constant  within  the 
Wigner-Seitz  cell.  The  e-density  in  RM  differs  from 
zero  around  the  boundary  of  the  Wigner-Seitz  cell. 
According  to  the  results  of  numerical  simulations,  the 
behavior  of  e-density  and  especially  the  behavior  of  a 
single-electron  potential  (e-potential)  are  quite 
unusual.  The  e-potential  displays  a clearly  pronounced 
minimum  wherever  the  e-density  reaches  clearly  pro- 
nounced maxima!  This  is  due  to  the  fact  that  the 
exchange  interaction  brings  electrons  to  the  area  of 
increased  e-density.  A similar  effect  has  been  known 
for  a long  time  in  the  theory  of  metals.  This  effect  is 
described  by  the  Slater  exchange  potential  Vexc(r)  = 
-a[p(r)]l/3,  where  p(r)  is  the  e-density  at  the  point  r 
and  a = 2 eV  is  a constant.  The  contribution  of  this 
effect  in  usual  metals  is  normally  several  percent.  In  the 
case  of  RM  properties,  this  effect  plays  a crucial  role. 

This  important  conclusion  of  CES  theory  is  a phys- 
ical manifestation  of  the  compensation  of  e-e  interac- 
tion due  to  the  identity  of  electrons  as  Fermi  particles 
whose  quantum  state  in  RM  is  described  by  an  antisym- 
metric wave  function. 

More  detailed  calculations  show  that  conductivity 
electrons  are  accumulated  around  the  minima  of  the 
e-potential  in  the  form  of  a degenerate  electron  gas,  giv- 


ing rise  to  a potential  barrier  for  the  penetration  of  con- 
ductivity electrons  to  the  center  of  the  Wigner-Seitz 
cell.  Thus,  the  exchange  interaction,  on  the  one  hand, 
assists  the  stabilization  of  the  system,  on  the  other 
hand,  may  considerably  increase  the  lifetime  of  the 
metastable  RM  state  up  to  macroscopic  time  intervals 
on  the  order  of  several  minutes.  The  latter  effect  was 
observed  experimentally. 

The  experimental  results  have  been  obtained  using 
new  techniques  to  create  large  densities  of  highly 
excited  Rydberg  species.  Alkaline  atom  Rydberg  spe- 
cies are  generally  simplest  to  form  because  of  the  low 
ionization  energy  of  alkaline  atoms.  The  new  tech- 
niques employ  diffusion  of  alkaline  ions  from  the  bulk 
of  nonmetallic  materials.  For  example,  in  the  case  of 
graphite  surfaces,  it  was  shown  that  the  emission  of  Cs+ 
ions  gives  rise  to  a nearly  resonant  process  that  forms 
Rydberg  species  from  the  emitted  ions  and  thermal 
electrons.  By  increasing  the  pressure  of  Cs  vapor  in 
contact  with  such  surfaces  and  by  ensuring  that  diffu- 
sion takes  place  in  the  material,  high  densities  of  Cs  and 
K Rydberg  species  have  been  achieved.  Small  particles 
(clusters)  and  macroscopic  amounts  of  RM  can  be  pro- 
duced by  using  Cs  vapor.  Alkaline-doped  metal-oxide 
surfaces  (e.g.,  promoted  catalyst  surfaces)  have  been 
also  shown  to  provide  high  densities  of  K Rydberg  spe- 
cies. The  behavior  observed  for  RM  produced  by  such 
methods  qualitatively  agrees  with  theoretical  predic- 
tions: RM  is,  for  example,  transparent  to  the  visible 
light  and  emits  much  less  light  than  an  ordinary  plasma 
would  when  high  currents  are  carried.  The  observed 
matter  is  very  energetic,  which  is  easy  to  observe 
through  explosions  of  the  matter  under  simultaneous 
emission  of  the  visible  light  or  charged  particles  for 
both  small  particles  of  RM  and  large  layers  of  RM  on 
surfaces.  A more  detailed  discussion  can  be  found  in 
[10,  13]. 

Since  the  first  experimental  observations  of  RM  in  a 
cesium  metal  plasma,  spectroscopic  studies  of  this 
highly  excited  form  of  condensed  matter  have  been  per- 
formed [32].  Electronic  Raman  anti-Stokes  and  mixed 
Stokes  processes  have  been  detected  in  the  range  of 
290-920  nm  upon  nonresonant  excitation  at  564  nm 
using  a nanosecond  pulsed  dye  laser.  No  fluorescence 
was  observed  in  this  first  Raman  study  of  the  RM  form 
of  Cs.  Stokes  transitions  have  been  mixed  in  one  more 
step  with  laser  radiation  and  go  from  low  circular  Cs+ 
states  to  RM  states.  Thus,  it  was  concluded  that  Raman 
spectroscopy  in  general  should  be  the  most  useful 
method  to  study  RM  due  to  the  extremely  large  polar- 
izability of  highly  excited  Rydberg  states,  forming  the 
RM  phase  by  condensation. 

It  was  shown  recently  [33]  that  much  less  powerful 
radiation  of  cw  IR  lasers  also  interacts  with  RM 
through  stimulated  Raman  processes,  revealing  novel 
features  of  RM  and  opening  more  possibilities  for  sys- 
tematic studies  of  various  density  regimes  of  RM  exist- 
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ence.  It  was  found  that  the  wave-number  blue  shifts  of 
single-mode  laser  radiation  in  low-density  RM  are  due 
to  the  stimulated  Raman  process.  In  this  case,  RM  is 
mainly  produced  from  Rydberg  states  of  K atoms. 
These  experiments  were  performed  with  milliwatt  cw 
lead  salt  tunable  diode  lasers.  Blue  shifting  was 
observed  as  shifts  of  Fabry-Perot  etalon  fringes  and 
was  on  the  order  of  0.05  cm-1.  This  blue  shifting  is 
attributed  to  a complete  conversion  to  anti-Stokes  radi- 
ation through  higher  order  stimulated  Raman  shifting 
by  extremely  small  quanta  from  the  near  continuum  of 
energy  levels  in  the  metallic  Rydberg  matter  phase.  The 
high  efficiency  of  the  stimulated  Raman  process  is 
believed  to  be  partially  due  to  the  weak  coupling  with 
other  degrees  of  freedom. 

Recent  experiments  have  also  demonstrated  stimu- 
lated anti-Stokes  Raman  scattering  from  an  RM  surface 
boundary  layer  on  the  sample  of  the  same  type  as  the 
one  used  in  [33].  These  experiments  were  performed 
with  a low-intensity  He-Ne  laser  with  a spectrometer 
used  to  analyze  the  light  coming  from  the  sample  anti- 
parallel to  the  laser  beam.  The  anti-Stokes  lines  from  K 
atoms,  as  well  as  from  H2  and  other  small  molecules 
were  also  observed  along  with  RM  transition  bands 
[34],  The  corresponding  Stokes  features  were  not 
observed,  which  is  a direct  indication  of  the  inversion 
of  the  system,  in  agreement  with  assumptions  and  dem- 
onstrations of  [33]. 

Rydberg  matter,  i.e.,  a condensed  phase  of  excited 
atoms  and  molecules,  was  predicted  more  than  20  years 
ago.  Comprehensive  quantum-mechanical  calculations 
have  been  performed  to  predict  different  properties  and 
the  lifetime  of  RM.  These  calculations  have  shown  that 
RM  is  a long-lived  excited  state  of  matter.  According  to 
theoretical  predictions,  the  RM  lifetime  may  be  as  long 
as  10  min  or  even  several  hours  and  days.  These  predic- 
tions are  consistent  with  long  radiative  lifetimes  of  even 
isolated  circular  Rydberg  atoms.  (The  lifetime  aver- 
aged over  the  angular  momentum  quantum  number  is 

0.1 8 s for  n = 40  and  17  s for  n = 100  [32].)  De-excita- 
tion was  found  to  occur  mainly  due  to  the  Auger  pro- 
cesses [10]  involving  two  electrons,  which  simulta- 
neously change  their  energies  and  orbital  momenta.  It 
was  demonstrated  in  [33]  that  RM  has  a rather  long 
relaxation  time  for  electronic  excitation  with  a lower 
limit  of  10  ms. 

Rydberg  matter  may  be  involved  in  many  natural 
phenomena  both  on  Earth  and  in  the  interstellar  space. 
The  hypothesis  that  RM  may  be  very  common  in  inter- 
stellar space  has  been  accentuated  by  studies  showing  a 
good  agreement  with  the  so-called  unidentified  infrared 
(UIR)  bands,  which  dominate  the  emission  spectrum 
from  interstellar  space  [35].  Rydberg  matter  is  on  the 
borderline  between  the  classical  and  quantum-mechan- 
ical descriptions  of  matter  and  is,  therefore,  of  consid- 
erable fundamental  importance.  In  this  context,  Holm- 
lid  [36]  proved  recently  that  all  of  the  highly  excited 
electrons  have  to  more  in  one  plane.  Thus,  RM,  at  least 


in  its  highly  excited  states,  has  a form  of  monolayers  of 
atoms.  Because  of  the  low  surface  work  function  of  RM 
(according  to  theoretical  predictions  and  experimental 
results  [13],  its  values  fall  within  the  range  of 
0.2-0.8  eV),  many  applications  become  possible  in  the 
area  of  power  electronics,  energy  conversion,  optoelec- 
tronics, and,  perhaps,  quantum  computing. 
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